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Abstract: We prove Gibbs distribution of two-state spin 
systems(also known as binary Markov random fields) without 
hard constrains on a tree exhibits strong spatial mixing(also 
known as strong correlation decay), under the assumption that, 
for arbitrary 'external field', the absolute value of 'inverse 
temperature' is small, or the 'external field' is uniformly 
large or small. The first condition on 'inverse temperature' 
is tight if the distribution is restricted to ferromagnetic or 
antiferromagnetic Ising models. 

Thanks to Weitz's self-avoiding tree, we extends the result 
for sparse on average graphs, which generahzes part of the 
recent work of Mossel and Sly[15], who proved the strong 
spatial mixing property for ferromagnetic Ising model. Our 
proof yields a different approach, carefully exploiting the 
monotonicity of local recursion. To our best knowledge, the 
second condition of 'external field' for strong spatial mixing in 
this paper is first considered and stated in term of 'maximum 
average degree' and 'interaction energy'. As an application, 
we present an FPTAS for partition functions of two-state spin 
models without hard constrains under the above assumptions 
in a general family of graphs including interesting bounded 
degree graphs. 

Keywords: Strong Spatial Mixing; Self-Avoiding Trees; 
Two-State Spin Systems; Ising Models; FPTAS; Partition 
Function 

I. Introduction 

Counting problem has played an important role in theoretic 
computer science since Valiant[19] introduced #P-Complete 
conception and proved many enumeration problems are com- 
putationally intractable. The most successful and powerful 
existing method for counting problem is due to Markov 
Chain method, which has been successfully used to pro- 
vide a fully polynomial randomized approximation schemes 
(FPRAS)(which approximates the real value within a factor 
of e in polynomial time of the input and with the 
probability > 3/4) for convex bodies[3] and the number of 
perfect matchings on bipartite graphs [9]. Since many counting 
problems such as the number of matchings, independent sets. 



circuits[14] etc. can be viewed as special cases of computing 
partition functions associated with Gibbs measures in statis- 
tical physics. Hence studying the computation of partition 
function is a natural extension of counting problems. 

Self-reducing [10] or conditional probability method is a 
well known method to compute partition functions if the 
marginal probability of a vertex can be efficiently approxi- 
mated. Gibbs sampling also known as Glauber dynamics is 
a popular used method to approximate marginal probability. 
This is a Markov Chain approach locally updating the chain 
according to conditional Gibbs measure. Hence studying the 
convergence rate(also known as mixing time) of Glauber 
dynamics becomes a major research direction. Recently the 
problem whether the Glauber dynamics converges 'fast' (in a 
polynomial time of the input and logarithm of reciprocal of 
sampling error) deeply related to whether a phase transition 
takes place in statistical model has been extensively studied, 
see [16] for hard core model(also known as independent set 
model) and [15][6] for ferromagnetic Ising model. Another 
approach to approximate marginal probability comes from the 
property of the structure of Gibbs measures on various graphs. 
This method utilizes local recursion and leads to deterministic 
approximation schemes rather than random approximation 
schemes of Markov Chain method. Our paper focuses on this 
recursive approach. 

The recursive approach for counting problems is introduced 
by Weitz[21] and Bandyopadhyay, Gamarnik [1] for counting 
the number of independent sets and colorings. The key of this 
method is to establish the strong spatial mixing property 
also known as strong correlation decay on certain defined 
rooted trees, which means the marginal probability of the 
root is asymptotically independent of the configuration on the 
leaves far below. Usually the exponential decay with the dis- 
tance implies a deterministic polynomial time approximating 
algorithm for marginal probability of the root. In [21], Weitz 
establishes the equivalence between the marginal probability 
of a vertex in a general graph G and that of the root of a 
tree named self -avoiding tree associated with G for two- 
state spin systems and shows the correlations on any graph 



decay at least as fast as its corresponding self-avoiding tree. He 
also proves the strong correlation decay for hard-core model 
on bounded degree trees. Later Gamarnik et.al.[5] and Bayati 
et.al.[2] bypass the construction of a self-avoiding tree, by 
instead creating a certain computation tree and estabUshing 
the strong correlation decay on the corresponding computation 
tree for Ust coloring and matchings problems. An interesting 
relation between self-avoiding tree and computation tree is that 
they share the same recursive formula for hard-core model. 
Considering the motivation of construction of the self-avoiding 
tree, Jung and Shah[8] and Nair and Tetali [17] generaUze 
Weitz's work for certain Markov random field models, and Lu 
et.al.[12] for TP decoding problem. Mossel and Sly[15] show 
ferromagnetic Ising model exhibits strong correlation decay on 
'sparse on average' graph under the tight assumption that the 
'inverse temperature' in term of 'maximum average degree' is 
small. 

In this paper, based on self-avoiding tree, we establish the 
strong spatial mixing for general two-state spin systems also 
know binary Markov random field without hard constrains on 
a graph that are sparse on average under certain assumptions. 
Our first condition is on the 'inverse temperature'. We show 
that there exits a value Jd in term of 'maximum average 
degree' d, if the absolute value of the 'inverse temperature' 
is smaller than Jj, for arbitrary 'external field', the Gibbs 
measure exhibits strong spatial mixing on a sparse on average 
graph. Since for (anti)ferromagnetic Ising model, strong spatial 
mixing on a finite regular tree implies uniqueness of Gibbs 
measures of infinite regular tree[4][20]. Jd in our setting 
is the critical point for uniqueness of Gibbs measures of 
infinite regular tree with degree of each vertex d, implying 
our condition is also necessary on trees. The condition is 
the same as that of Mossel and Sly[15] when ferromagnetic 
Ising model is the only focus. This makes part of their work 
in our framework. Our proof yields a different approach, 
and also avoids the argument between weak spatial mixing 
and strong spatial mixing employed in |21]. In fact our 
proof is based an inequahty similar to the one in [13] and 
carefully exploits monotonicity of the recursive formula. The 
recursive formula on trees is well known. Recently Pemantle 
and Peres [18] use it to present the exact capacity criteria 
that govern behavior at critical point of ferromagnetic Ising 
model on trees under various boundary conditions. Our second 
condition is for 'external field'. We prove for any 'inverse 
temperature' on a graph which are sparse on average Gibbs 
distribution exhibits strong spatial mixing when the 'external 
field' is uniformly larger than -B(ci, amax, l) or smaller than 
—B{d, — amin, 7), where d is 'maximum average degree' and 
ctmin, ctmax, 7 are parameters of the system. To our best 
knowledge, this condition on 'external field' is first considered 
for strong spatial mixing. The technique employed in the proof 
is Lipchitz approach which has been used in [1][2][5]. The 
novelty here is that we propose a 'path' characterization of 
this method, allowing us to give the 'external field' condition 
in term of 'maximum average degree' rather than maximum 
degree. Some notations of the 'sparse on average' graphs have 



appeared in [15]. These are graphs where the sum degrees 
along each self-avoiding path(a path with distinct vertices) 
with length O(logn) is O(logn). 

As an application of our results, we present a fully poly- 
nomial time approximation schemes(FPTAS)(which approxi- 
mates the real value within a factor of e in polynomial time of 
the input and e~^) for partition functions of two- state spin 
systems without hard constrains under our assumptions on 
the graph G = {V,E), where, for each vertex v G V of 
G, the number of total vertices of its associated self-avoiding 
tree Tsaw(v) with hight O(logn) is 0(n'^(^)). This includes 
bounded degree graph and especially Z'^ lattice more con- 
cerned in statistical physics. Jerrum and Sinclair [7] provided 
an FPRAS for partition function of ferromangetic Ising model 
for any graph with any positive 'inverse temperature' and 
identical external field for all the vertices. Their results do 
not include the case where different vertices have different 
external field, and are not applied to antiferromagnetic Ising 
model where the 'inverse temperature' is negative either. 

The remainder of the paper has the following structure. In 
Section II, we present some preliminary definitions and main 
results. We go on to prove the theorems in Section 111. Section 
IV is devoted to propose an FPTAS for the partition functions 
under our conditions. Further work and conclusion are given 
in Section IV. 

II. Preliminaries and Main Results 

A. Two-State Spin Systems 

In the two-state spin systems on a finite graph G = {V, E) 
with vertices V = {1,2, and edge set E, a 
configuration consists of an assignment cr = (cti) of 
D, = {±1} values, or "spins", to each vertex(or"sites") of 
V. Each vertex i G V is associated with a random variable 
Xi with range ±1. We often refer to the spin values ±1 
as (+) and (— ). The probability of finding the system in 
configuration cr € fi" is given by the joint distribution of 
n dimensional random vector X = {Xi,X2, - ■ ■ ,X„}(also 
known as the Gibbs distribution with the nearest neighbor 
interaction) 

Pg{X = a) = exp( ^ /3y (fji, Uj) + ^ hiipi)). 

Here Z{G) is called partition function of the system and a 
normalized factor such that X^cren" Pg{X = a) = 1, and hi 
and /3ij are defined as a function from fl and fi^ to i?U{±oo} 
respectively. We use notation (3ij{a, h) = (3ji{b, a). We say the 
system has hard constraints if there exit an edge G E 
or a vertex k, and an assignment ai, aj or ctj. such that 
Pij{<^i\(^j) = oo OT tikijJk) = 00 (e.g. hard-core model is one 
of the systems with hard constrains where +) = —oo, 

-) = = -) = and /7,:(-) = 0). In 

this paper we focus to the systems without hard constrains. 
We call the function 'interaction energy' and hi 'appUed 
field' . If /3ij((Ti, (jj) ~ Jijaiaj and hi — Bicri for all the edge 
G E and vertex i G V, where and Bi are constant 



numbers varying with edges or vertices, the system is called 
Ising model. Further, if is uniformly (negative )positive 
for all e E, the system is called (anti)ferromagnetic 

Ising model. Jij and Bi are called inverse temperature 
and external field respectively. To match the notation of 
Ising model, set J„- = 

and Bi — '^'M^'^'^^') for all edges and vertices , in this paper 
we call Jij and Bi are 'inverse temperature' and 'external 
field' of general two-state spin systems without hard constrains 
(denoted by TSSHC for abbreviation)respectively. For any 
A C V, aA denotes the set {ai,i E A}. With a little abuse 
of notations, cta also denotes the configuration that i is fixed 
(Ti, Vi E A. Let Z{G, $) denote the partition function under 
the condition $, e.g. Z{G, Xi — +) represent the partition 
function under the condition the vertex 1 is fixed +. 

B. Definitions and Notations 

Definition 2.1 (Self-Avoiding Tree) Consider a graph 
G — (V, E) and a vertex v E V in G. Given any order of all 
the vertices in G. There is associated partial order on E of 
the order on V defined as > (k.l) iff {i,j), {k,l) share 
a common vertex and i + j > k + I. The self-avoiding tree 
Tsaw{v){G){foi simplicity denoted by T5a^(^)) corresponding 
to the vertex v is the tree of self-avoiding walks originating 
at v except that the vertices closing a cycle are also included 
in the tree and are fixed to be either + or — . Specifically, 
the vertex of the Tgaw{v) closing a cycle is fixed -I- if the 
edge ending the cycle is larger than the edge starting the 
cycle and — otherwise. Given any configuration cta of G, 
A C 1^, the self-avoiding tree is constructed the same as the 
above procedure except that the vertex which is a copy of the 
vertex i in A is fixed to the same spin ai as i and the subtree 
below it is not constructed(See Figure 1). Hence, for any 
configuration of G, A C V, we also use cta to denote the 
configuration of r^aii;(t>) obtained by imposing the condition 
corresponding to cta as above. For any E E and i E V 
of G, the 'interaction energy' function and 'applied field' 
function on all their copies of the induced system on Tsaw(v) 
by G are the same as Pij and hi respectively. 

We now provide the remarkable property of the self- 
avoiding tree, one of two main results of [21], which is one 
of the essential techniques of our proofs. 
Proposition 2.1 For two-state spin systems on G — (V, E), 
for any configuration a a, A C V and any vertex v EV, then 

PaiXy = +|aA) = -Pt,_(„,(A:. = +|cta). 



In order to generalize our result to more general families 
of graphs , which are sparse on average, we need some 
definitions and notation of these graphs. 

Definition 2.2 Let \A\ denote the cardinality of the 
set A. The length of a path is the number of edges it 
contains. The distance of two vertices in a graph is the 
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Fig. 1. The graph with one vertex assigned + (Right)and its 
corresponding self-avoiding tree rg„^(i)(Left) 

length of shortest path connecting these two vertices. A 
path vi,V2,--- is called a self-avoiding path if for all 
* 7^ J ' ''^i 7^ ^'j- In a graph G = {V,E), let d{u,v) denote 
the distance between u and v, u,v E V. The distance 
between a vertex v E V and a subset A C V is defined by 
d{v, A) =min{d('i;, u) : u E A}. The set of vertices within 
distance ^ of w is denoted by V{G, v, I) ^ {u : d{v, u) < I}. 
Similarly, the set of vertices with distance Z of u is denoted by 
S{G, V, I) = {u : d{v, u) = I}. We call a vertex at the height 
t of a rooted tree if the distance between it and the root is 
t. Let Sy denote the degree of v in G. The maximal path 
density m is defined by m(G, v, I) = max ^ S^, where the 

maximum is taken over all self-avoiding paths T starting at v 
with length at most I. The maximum average path degree 
S{G,vJ) is defined by SiG,vJ) = im{G,v,l) - Sy)/l,l > 1. 
The maximum average degree of G is defined by 
A{G,l) — maxy^v S{G,v,l). Roughly speaking, in this 
paper, a family of graphs Q is sparse on average if there exits 
a constant number a and d such that A(G, alogn) < d for 
any G eQ. 

Some properties of the above definitions are useful in our 
proof, we present them. Most of proofs are simply obtained 
by induction and can be found in [15]. 

Proposition 2.2 Let j, I denote positive natural numbers, 
then 

m{G, v,jl) < j max{m(G, u, I) — dy} + Sy. 

MSG 



Proposition 2.3 Let I be natural numbers, then 

|5(T,,^,(,),i;,Z + 1)1 < 5y[5{G,v,l) - l)K 



Proposition 2.4 Let j, I be natural numbers, then 

\V{Tsaw{v),v,il)\ < {max\V{Tsaw{u),u,l)\y . 



Definition 2.3 ((Exponential) Strong Spatial Mixing) Let 
G = {V, E) be a graph with n vertices. The Gibbs distribution 
of two-state spin systems on G exhibits strong spatial mixing 
iff for any vertex v gV, subset A cV, any two configurations 
CTA and rjA on A, denote Q = {v G A : ^ r]v} and 

t = d{v,e), 

\PG{Xy = +|(7a) - Pg{X, = +|77a)| < fit), 

where f{t) goes to zero if t goes to infinity and is called decay 
function. 

For the purpose of our settings, we present a weak form 
of exponential strong spatial mixing. We say the distribution 

exhibits exponential strong spatial mixing if there exits positive 
numbers a, b, c independent of n such that f{t) = bexp{—ct) 
when t = kalogn, k = 1,2, - ■ ■ . 

Remark: In the above definition of (exponential) strong 
spatial mixing, PciXy = +|crA) and Pa^Xy = +|?7a) can be 
replaced by \og{PG{Xy = and \og{PG{Xy = +\r]A)) 

respectively if d{v, A) is large than a constant number, due 
to the inequality 2x <log(l + x) < x when |a;| < .5, and we 
call it the logarithmic form exponential strong spatial mixing. ) 

Definition 2.4 (FPTAS) An approximation algorithm is 
called a fully polynomial time approximation scheme(FPTAS) 
iff for any e > 0, it takes a polynomial time of input and 
to output a value M satisfying 



where M is the real value. 

Remark: In the above definition 1 — e and 1 + e can be 
replaced by and e^. 

C. Main Results 

For simplicity , We use the following notations. Consider 
a two-state spin systems with hard constrains(TSSHC) on 
a graph G = {V,E) with ri vertices V = {1,2, ••• ,n} 
and edge set E. Let J = max(jj)g£; \ Jij\, 

Brain = minjgy -Bmax = max,gyi?i, ttniax = 

max {A,(-,-) - - 
amin = min - l3ij {+,-), Pi - 

fi..(j_ J-W TTiQY/ N r \bijCij — aijdij\ \bijCij — aijdij \ i 

7 = max(jj)£B{7y}' where Jy = 

o _ hi(+)-hi(-) 

are 'inverse temperature' and 'external field' respectively, 

and fly cxp(/3ij (+,+)), = exp(/3jj(-F, -)), 

Cij = exp(/3y (-, +)), dij = exp{Pij{-, -)). 

Theorem 2.1 Let G = {V, E) he a graph with vertices V — 
{1, 2, • • • , n}, edges set E and TSSHC on it. If there exit two 



positive numbers a > and d > such that A(G, alogn) < 
d, and when 

{d- l)tanhJ < 1 

or equivalently J < J4 = ^^og{-^^), then the Gibbs 
distribution of TSSHC exhibits logarithmic form exponential 

strong spatial mixing for arbitrary 'external field', specifically, 
for any i € V, any two configurations a\ and rj\ on A, denote 
@ = {j G A : aj r]j} and t = d{i,Q) — kalogn + 1, 
k = l,2,---, 

I l0g{PG{Xi = +|c7a)) - log{PG{Xi = +|r?A))| < fit), 
where f{t) = 4J^i((d- l)tanh J)*-i. 

Remark: If the graph is bounded with the maximum degree 
D, then d can be replaced by D while for any a > 0, 
and Jij is the 'inverse temperature' in (anti)ferromagnetic 
Ising model, then theorem 2.1 still holds and Jd is the 
critical point for uniqueness of Gibbs measures on a infinite 
tree with maximum degree _D[13]. Note the decay function 
is slight different from the definition since 5i may be 
O(logn), however, in this case we can choose k large enough 
independent of n such that f{t) = e"*"* when n is large, 
where & is a positive number independent of n, then replace a 
by ka as required. In fact in the apphcation of the algorithm, 
this is not important. 

Theorem 2.2 Let G = {V, E) be a graph with vertices V = 
{1, 2, • • • , n}, edges set E and TSSHC on it. If there exit two 
positive numbers a > and d > such that A(G, alogn) < 
d, and [d — 1) tanh J > 1, and when 

7) or Bmax < -amin,7) 

Where B{d,a,j) = + log( ^^^^+/^^" -^ibf ), the 

Gibbs distribution of TSSHC exhibits exponential strong spa- 
tial mixing, specifically, for any i € V, any two configurations 
(TA ond ?7a on A, denote = {j G A : Uj ^ 77^} and 
t = d{i,e) = ka\ogn+l. A: = 1,2, • • • , 

|Pg(X, = +|aA) - PGiX^ = +|r,A)| < fit), 

where fit) = '^'^ ( (rf-l)7exp(2B„i„-(d-l)a„^^) w-1 
fit) - 'Si7 ( (d-l)7exp(2B„ax-(rf-l)a„in) u-1 resDcctivelv 

Jyt) - 4 K (i+exp(2B„.ax-(d-i)a„i„))2 > respectively 

Remark: It's easy to check 7 > 4 tanh J, hence in theorem 
2.2, if id - 1) tanh J > 1, then 7(d - 1) - 4 > 0. As 
a corollary of Theorem 2.2, from its proof in section III, 
we know if the graph is bounded degree with maximum 
degree is d, the condition for 'external field' can be relaxed to 
Bi > Bid, amax, 7) or Bi < -Bid, -a^in, 7) for anyiGV, 
which does not require that 'external field' is uniformly large 
or uniformly small. 

Theorem 2.3 Let G = (V, E) he a graph with n vertices 
V = {1, 2, • • • , n}, edges set E and TSSHC on it. If there 



exit two positive numbers a > and d> Q such that for any 

i e V 

t/(r,„^(,),i,alogn)<(ci-iri°8", 

where \V{Tsawii),i,l)\ = {j e Tsaw(i) ■ d{i,j) < I}, 
then when J < or J > J^, -Bmin > B{d, a^a.x,j) or 
-Bmax < —B{d, — amin, 7), there exits an FPTASfor partition 
function of TSSHC on G. 

III. Proofs 

We now proceed to prove the theorems. One of the 
technical lemmas for the theorem 2.1 is an inequaUty similar 
to [13]. We present it now. 

Lemma 3.1 Let a, b, a, d, x, y be positive numbers and 



9{x) = 



ad—^/hc 



/ ad+^/bc 



, then 



max( 



9{x) g{y) 
9{y) ' 9{x) 



)<(max(-,^))*. 
y X 



Proof: Case \.ad> he. Since g{x) = 



ad— be 



cx+d c c{cx+d) 

is an increasing function, w.l.o.g. suppose x > y and let x = 
zy, where z > 1, then 



log( 



9M) = r5!!!(5i^ 

9{y) Ji da 



ay 



da 

cy 



-I 



)da 

1 {aay + b) (cay + d) ^'^ 

(ad — hc)y 



aay + b cay + d 
{ad — bc)y 



< 



1 {■•/acay — \/bd)'^ + {Vbc + \fad)'^ay 
{ad — bc)y 



da 



(v 6c + V adYay 



^/ad — \/bc , 
da = — logz. 



ad - 



be 



Hence, 



,q{x) q{v) . q{x) ,x 

9{y) 9{x) g{y) y 



(max(-,^))*, 

y X 



where t = 

V aa+v oc 

Case 2. ad < be. Similar to the first case, g{x) is a decreasing 
function, let h{x) = l/g{x), then h{x) is an increasing 
function, w.l.o.g. suppose x > y, repeat the process of Case 1 
for gg, then 



%) - ^y' 



Hence, 



max( 



where t 



9ix) g{y) ~^ ^ g{y)_ ^ h{x)_ ^ ^xy 

9{y) ' 9{x) g{x) h{y) ~ y 
□ 



(max(-,^))*, 

y X 



i ad 



/ad+Vbc' 



Lemma 3.2 Let T = {V,E) be a tree rooted at with 
vertices V = {0, 1,2,--- edge set E and TSSHC 



on it. Suppose some vertices are fixed. Let Tk and Ti be 
two subtrees of T including vertex k and I respectively by 
removing an edge {k,l) where d{k,0) < d{l,0). The fixed 
vertices remain fixed on Tk and Ti. Then the probability of 
Xq = + onT equals the probability of Xq = + on the subtree 
Tk except changing the 'external field' hk to certain value ft^.. 

Proof: Let CIti denote the configuration spaces, Ei and Vi 
the edge set and vertices on T(. Setting 

/ife(o-fc) = hk{ak)+ 

log( e <-^'^^' ) 

completes the proof. □ 

With Lemma 3.1 and Lemma 3.2, we now proceed to 
prove (exponential) strong spatial mixing property on trees. 

Theorem 3.1 Let T be a tree rooted at with vertices 
y = {0, 1, 2, • • • ,n}, edge set E and TSSHC on it Let AcV 

, and 77a be any two configurations on A. Let Q = {i : 
Ci i- r]i,i ^ A}, t = d(0,e) and s = \S{T,0,t)\ = \{i : 
d{0,i) =t,i€ T}\. Then 



max( 



Pt{Xo = +|Ca) Pt{Xo = 
Pt{Xo = +\Va)' Pt{Xo = 



) < exp(4Js(tanh J)*~^) 



Proof: For any i G V, let Ti denote the subtree with i 
as its root and Z{i) be the TSSHC induced on by T. 
Noting To is T. To prove the theorem, it's convenient to deal 
with the ratio pl[xlZt\'ll] rather than Pt{Xo = +|Ca) itself. 



Denote R, 



where ^Ai is the condition by 



imposing the configuration ^a on Ti, and note a simple relation 
if XUX2 e (0,1), then |i > l iff ^ > further 

niax{^,2^} < max{24T^,^4T^}- Hence replace 



xi and X2 by Pt {Xo 
only to show 



+ |Ca) and Pt{Xo = +|?7a), we need 



^^i^il' —77) ^ exp(4Js(tanh J)*"^). 

-"'0 Rq 



(1) 



Theorem 3.1 follows by max( 



Pt{X„ = + \Ca) 



Pt(Xo = + |j)a) ' Pt(Xo=+|Ca) 



Pr(Xo=+|r)A) ^ <• 



max(^,^). 
We go on to prove (1) by induction on t. Before we doing 
this, some trivial cases need to be clarified. We are interested 
in the case t>l and is unfixed. Let Tki denote the unique 
self-avoiding path from fc to Z on T. If i is a leave on T and 
d{0,i) < t, where t = d{0,e). Define U = {j G V : j G 
Toi,3k e S{T,0,t),s.t.j e Tofe}. Note U 9 since G U. 
Let ji E U such that d{i,ji) = d{i,U). By lemma 3.2, we 
can remove the subtree bellow ji and change external field 
hj^ at ji to hj. without changing the probability of Xq = +. 
More importantly, this process removes at least one leave at 
the hight < t, and does not remove any vertex at the hight > t. 



Thus, w.l.o.g. suppose T is a tree rooted at where any leave 
on it at the height > t. Let Oi, O2, • • • > Og be the neighbors 
connected to 0. A trivial calculation then gives that 



^Ca ^ Z{To,Xo = +, Ca) 
° Z(To,Xo = -,Ca) 



where the last equation follows by '}2'-=i — ^- This 
completes the proof of Theorem 3.1. □ 

With Theorem 3.1 and self-avoiding tree, it's enough to 
prove Theorem 2.1. 



ID(/3oOi(-,o-Oj)+ J2 Pkl{o-km)+ J2 hk(<Tk)) 

creQo 

/3oOi(+,<^Oi)+ J2 0kl{crk,(7i)+ J2 hk((7k) 

g (k,l)€Ti kSTf 



n 



„2Br 



= g2Bo 



n 

i=l 

n 



-,Ca,) 



,^1 CiZ{To,,X^ = +, CaJ + diZ{To-,Xi = -, CaJ 



dCa 



(2) 



where = ^°(+)-^°(-) , = e'3oo.(+.+), 6^ = e'3oo.(+.-), 
a = e^°"i'^~'+\ di = e''™i(~'~) . Now we check the base 
case t = 1 where Rq^,R^'^ G [0, +00], by the monotonicity 



of 



and 



Ca p^a 1 



max(-; — , — -) < 



•0 



6^Cj (lid{ 



Proof of Theorem 2.1: Due to Proposition 2.1, the 

only thing left is to verify \S{Tgaw{i)iht)\ = Si{d — l)*"-"- 
when t = d{i,Q) = ka log n + 1, k = 1,2, under the 
condition that there exit two positive numbers a > and 
(i > such that A(G', alogn) < d. By Proposition 2.2, 
we know m{G,i,ka\ogn) < ka\ognA{G, alogn) + Si < 
kalognd + Si, hence, S{G, i, kalogn) < d follows from the 
definition S{G,i,l) = {■m{G,i,l) — Si) /I. By Proposition 
2.3, it's sufficient to show |S'(Tsoj(,(j), z, fca logn + 1)| < 
5^{5{G,^, kalogn) - l)fe«i°g" < Si{d - l)'=«i°s" = 
5i{d — This is exactly what we need. □ 

Next we will proceed to prove Theorem 2.2, we still use the 
recursive formula but with another form. The technique used is 
a well known method, Lipchitz approach. A 'path' version of it 
will be presented, which allow us to bound the 'external field' 
with maximum average degree. Before presenting it, we need 
some notation for simplicity. Let T = (y, E) be a tree rooted 
at with vertices 0, 1,2,--- ,n, edge set E and TSSHG 
on it. For each edge G E, recall the notation in main 

results, flij = e^-^(+'+\ bij = e'^'^-(+'-\ = e'^'i(-'+\ 



and di^j 
Define 



fij{x) 



g/3«(- Let Mij = a 



d - 



K.i 



MijX ■ 



di'i 



NiiX + hi 



hij{x) 



{MijX + dij){NijX + bij) 



Hence t = 1, (1) holds. Assume by induction that (1) holds for 
t — 1, and we will show it holds for t. Let Si = |5'(ro; ,Oi,t — 
1)\, i = 1,2, - ■ ■ ,q, stiU using above recursive procedure, then 







i=l 
1 







) 



< 



^Ca JJ'?A . ,/a~d~-^ 

^ ^»?A ' pCa 



i=l 

q 



— 11 ^ R'JA ' nCA ' 



i=l 



where the second inequality comes from the Lemma 3.1. 

^Ca j^^A 

According to the hypothesis of induction max(-^^, -^) < 



exp(4Jsi(tanh J)* ^), it's sufficient to show 

.Ri"- R?^ 



^0; 



max(-2^, -2-) < exp(4Jsi(tanh J)*"!) 
J^o Ri 



i=l 



exp(4Js(tanh J)*"^) 



Recall 

Q!max = max {/?„■(-, -) - - 

(+,+)}. 

Q!mm = ^min - Pi j {+,-), fiij{-, +) - 



7 = max(j j')g£{7ij}. For any i E V, let Ti denote the 
subtree with i as its root and Z{i) be the TSSHC induced 
on T, by T. Recall Bi = ^'(+)-^'(-) is the 'external field'. 



-2Bi 



denote A, 



path from i to j on T. 



and let be the unique self-avoiding 



Lemma 3.3 For any {i,j) G E, max \hij{x)\ < -jij 

a;£[0,l] 

Proof: The proof is technique and left to the appendix. 



With the above notations, we present a 'path' version 
Lipchitz approach. 

Lemma 3.4 Let A C V , Ca ond rj^ be any two configu- 
rations on A. Let Q = {i : Q ^ r]i,i € A}, t = d(0, Q) and 



S{T, 0, t) = {i: d{0, i)=t,i€ T}. Then j with d{0, j) < t - 1 or there is a vertex i on T with 

\r> rv oc\^ ii M d(0,i) < t — 1 being fixed, the contribution of the path to 

\Pt{Xo = + Ca) - Pt{Xo = +\Va)\ • Ca Va Ca Va n 

the summation is zero since p> — p- = Pj — p- = 0. 

< 7* 5^ 9i{zi){i — gi{zi)) Hence the remaining path with length t is in the set 

kes{Tfi,t) lero.i^k {Tok ■ k e ^(r.O,^)}. This completes the proof of lemma 

where gi{xi) = (1 + Aj H fujixn^))^^ and Xi is a '-' 

(i,ij)6Ti 

vector with elements xu- G [0, € V, and Zi are constant In order to prove the Theorem 2.2, we need the following 

vectors with elements in [0, 1]. lemma. 

Proof: For any i in T, let p^^'' = Ft, (X,; = +|Ca,) Lemma 3.5 Let Xi>0, i = 1,2, - ■■ ,n. Then 

and R'i'^ = p^'|^'~^j^^'j , where is configuration by „ 

restriction of Ca on Tj. Then we have the following equality J^(|l _l_ > _|_ " 



pi'^ = Pt{Xo = +IC 



1 1 



A 



1 1 Proof: The proof is technique and left to the appendix. 



1+^0 n „„„ 1+Ao n 



(o,o,)6T ''oo,«of+*'oo,. (o.o^jeT ^oojPo,''+«'ooj With Lemma 3.4 and 3.5, it is sufficient to prove Theorem 

= go{xo) 

where xq = ' ' ' 'Po^ )• First, note for any x = of Theorem 2.2: Following the notation of Lemma 

(,ri, .7-2. ■ ■ ■ . :r„) and ?y = {jji^yl • • • , -7 = ^o, first order 3.4, let s = \S{T, 0, t)\, we have 
Taylor expansion at y gives that there exists a G [0, 1] such r/Ai ^ t TT 

feeS(T,0,t) ieFofci^fc 

5o(a:)-5o(t/) = V5o(y + e(:r-y))(x-yf, < V ^ max ^ [] " 

where {x — y)^ denotes the transportation of the vector ierofci^fc 

(a; - y). Calculate the we have ^7* tt / n/, ^ \x 

^ max [[ gMKl-gi{z,)). 

„ Ao 11 /oo,(xj)( j 

<^9o{x) ^ _ j=i By Lemma 3.5, for each Tojk, (0, Oj) e T, fc e ^(T, 0, t), 

dx- ^ 

(1 + Ao n/oo,(ari))^ n 9i{zi){l-gi{zi)) 

= -go{x)il - 9.(SC)){^^^1 - A, n 

_ TT eTi 

= go{x){l-goix)) .^^^ ^ . ^ . . ,-^r (1 + Ai 11 fikizin)? 



(i,ii)6T, 



Hence, there's exits € [0, 1] such that 

q where r^fc = ( H Aj f] /u, (^ii,))^/^*"^^- A simple 

^ ' calculation gives that e""'" < jij{x) < e""*'', for any S 
T. Hence, 
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<^±goizo){l - goizoM': "<! ^ e«-(^(^-°'*-i)-^). 

J=i Now we prove the exponential strong spatial mixing under 

where zo = p'i^ + 0o{pi^ - P^^) and the second inequality assumption of Theorem 2.2. Suppose T is a self-avoiding tree 

follows by Lemma 3.3. Now repeat the procedure for of G. (5(T, 0, 1) < A(G, t- 1) < d when t - fcalogn+ 1, 

bo^ - Potl j = 1,2, •••,9, it is easy to see that the k = 1,2,- ■■. If B^i^ > B{d,ama^,j), then 
summation is over all the self-avoiding paths emitting from lid — l)cxp(2i3min —ocraax{d— 1)) ^ 

the root 0. For each path T, if the end point of T is a leave (l + exp(2Smin - amax(<i - 1)))^ ^ 



Noting s < So{d- and ( 11 Aj)^/^*"^^ < e-^^"'" 
, now we can see 



Ca 



4 H^ + rjk)' 



^ (5o7 , 7((i - 1) exp(2.Bniin - ama^jd 
- 4 ^ (l + exp(2S^i„-a„,ax(d-l)))2 



1)))* 



(3) 



The similar case holds for Bmax < — -B(rf, — amin, 7)- This 
completes the proof. □ 

Remark: As we point out in section II that if the graph G is 
a bounded degree graph witii the maximum degree d, the con- 
dition in Theorem 2.2 can be relaxed to Bi > B{d, ctmax, 7) 
or Bi < — — amin, 7) for any i G V. The reason for 
this comes from the upper bound for gi{zi){l — gi{zi)) in 
the Lemma 3.4 since j{d — l)gi{zi){l — gi{zi)) < 1 for 
any i G To^k.i 7^ k where (0,0^) eT,ke S(T,0,t). We 
emphasize that one way to improve the condition by this 
method is to carefully analyze the bound of ,fij{x) for each 
iterative step according to the range of x since this will give 
better bound for gi{x). We do not optimize the parameter 
here and do not know whether dealing with the bound of 
fij{x) carefully makes tiie B{d, amax, 7) or -B{d, -a^in, 7) 
optimally approximate the critical point of 'external field' for 
uniqueness of Gibbs measures either if there does exit one(note 
that the critical points of 'external field' for ferromagnetic and 
antiferromagnetic Ising model are different on Cayley free, an 
infinite regular tree with the same degree for each vertex [4]). 

The proof of Theorem 2.3 will be shown in Section IV. 

IV. Approximating Partition Function 

In the proof of Theorem 2.1 and 2.2, the calculation of 
the marginal probabihty of the root yields a local recursive 
procedure. If we truncate the tree at height t, and then use 
the recursive method to compute the marginal probability at 
root, it is easy to see the complexity of this procedure is the 
number of vertices of truncated tree. We now present the 
algorithm based on the above procedure and self-avoiding 
tree. 

Let G = {V, E) be a graph with vertices F = {1, 2, • • • , n}, 
edge set E and TSSHG on it. Let $1 denote the whole state 
space ( which means Pg{X\ = -|-|$i) = Pg{Xi = +)), and 
= {Xi = +, 1 < i < j - 1}, 2 < j < n + 1. 



Algorithm for Partition Function Z{G) 

Input: G with the TSSHC, e > precision. 

Output: Z{G), the estimator of partition function Z{G). 

For j=1 :n 

compute pj, an estimator of conditional marginal 
probability pj = PciXj = through self-avoiding 

tree T^awU) truncated at a certain height tj under the 
condition $j such that (1-^) < < (1-^)- (The initial 
values of iteration at height tj are arbitrary nonnegative 



numbers, if we adopt the recursive formula in the proof of 
Lemma 3.4 where Pt{Xo = +\Ca) = 9o{xo)) 

Output:Z(G) = Z{G,^n+i) f[ Pi~'- 

i=l 

With the above algorithm, it is enough to prove Theorem 
2.3. 

Proof of Theorem 2.3: First we show under the as- 
sumption of the theorem, the Gibbs distribution exhibits 
exponential strong spatial mixing. Since(by Proposition 2.4) 

\V{Tsaw{i),i,ka\ogn)\ < \ max{V{Tsaw{]), j, ka\ogn))\'' < 

jev 

{d — l)'=''i°g" for any i e V and A; = 1, 2, • • • , we can obtain 
the trivial bound of the number of vertices at height ka log n, 
that is , \S{Tgaw(i),i, ka log n) < \V{T,aw{i),hka\ogn)\ < 
{d — log". Let t — kalogn from the proof of Theorem 
2.1 and 2.2(see Formula (1) and (2)), substituting {d — 1)* 
to s in (1) (2), we get the exponential strong spatial mixing 
of Theorem 2.3. Specifically, if J < J^, the decay function 
f{t) = 4:J{d — l){{d — l)tanh J)*~^ which corresponds to 
the logarithmic form exponential strong spatial mixing, and if 

> B{d 7) or Bmax <-B{d, -Q!,„i„, 7), 

the decay function has the same form as in Theorem 2.2 except 
replacing (5^ by rf— 1. In both cases, we suppose decay function 
f{t) = he^^^ where b, c are constant positive numbers inde- 
pendent oin,t = ka log n, fc = 1, 2, • • • . Through exponential 
decay property, it's sufficient to show the above algorithm 
provides an FPTAS for Z{G). Now we check the output Z{G) 
satisfying (1 



<(l + e). Since p, = ^±^. 



multiplying them gives Z{G) 



1 



Z(G,$„+i) n Hence, 

i=l 

we point out previously that the complexity of the algorithm at 
each step is 0{\V{T,^^^j),j,tj)\) = 0((rf- when tj = 
kalogn, k = 1,2,- ■ ■, we only need to set f{fj) < O(^) 
to promise (1 — ^) < < (■'-~^) which requires tj = 
0(log n + log(e^^)). Thus, the complexity of the algorithm is 

nO{(d~ l)0(logn+log(.-i))) ^ 0(„0(1) ^„(g-l)0(l))^ ^ijj^j^ 

completes the proof. □ 

V. Conclusion and Further Work 

We have shown that the Gibbs distribution of TSSHC on a 
'sparse on average' graph G = {V, E) with 'maximum average 
degree' d exhibits the (exponential) sttong spatial mixing when 
the absolute value of 'inverse temperature' | | < or the 
'external field' Bi is uniformly larger than ai„ax,7) or 
smaller than — — amin, 7), for any E E, i E V. 

Here Jd is the critical point for uniqueness of Gibbs measure 
on a infinite d regular tree of Ising model, implying the 
condition for inverse temperature is tight when restricting it 
on Ising model, B{d, a, 7) is constant with parameter d, a, 7. 
It is not difficult to apply our results to Erdo-Renyi random 
graph G{n,d/n), where each edge is chosen independently 
with probability d/n, since the average degree in G{n, d/n) 



is d{l — 0(1)) while it contains many vertices with degree 
log n/ log log n[15]. As an application of strong spatial mixing 
property, we present an FPTAS for partition functions on 
a little modified sparse graphs, which includes interesting 
bounded degree graph. 

For future work, we expect to improve the condition on 
'external field' . We have presented a way to improve it in the 
remark, however, we believe the essential improvement needs 
other method. Maybe the approach of analysis of the fixed 
point in[ll] works here. 
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Appendix 

Proof of Lemma 3.3: 

Since MijX + dij > and Nijx + bij > 0, Va; e [0, 1], 
we need only to show min w{x) = mm{aijCij,bijdij), 

x£[0,l] 

where w{x) = {MijX + dij){NijX + bij). The case 
MijNij = is trivial, so w.l.o.g. suppose MijNij ^ 0. 
Noting xi = — '''^ 2M^Jv'^'^ ™ extremum of w{x) on R. 
There are three cases needed to be discussed. 

Case 1. MijNij < 0, then 'w{x) reaches its minimum 

at boundary. Then min w{x) < min('u;(0), w(l)) = 
xe[o,i] 

m.m{aijCij,bijdij). 

Case 2. My > 0, Nij > 0, then xi < 0, 'w{x) is increasing 
on [0, 1], then min w{x) = w{0) = bijdij. 

xe[0,l] 

Case 3. My < 0, Nij < 0, then xi > 1, wix) is decreasing 

on [0, 1], hence min w{x) = w{l) = UijCij. □ 
xe[o,i] 

Proof of Lemma 3.5: 



n(i + Ao = i + E( E U^o 

fc=l n<i2<---<ifc j=l 



fe=l i=l 
n n 

=i+E(Cn(n^')*) 



fe=l i=l 



(1 



where = - The first inequality uses the arithmetic- 

geometric average inequaUty. 



